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The surface of a topological insulator is a closed two dimensional manifold. The surface states are 
described by the Dirac Hamiltonian in curved two dimensional spaces. For a slab-like sample with 
a magnetic field perpendicular to its top and bottom surfaces, there are chiral states delocalized on 
the four side faces. These "chiral sheets" carry both charge and spin currents. In strong magnetic 
fields, the quantized charge Hall effect {axy — (2n + l)e^ /h) will coexist with spin Hall effect. 



The theoretical prediction that 3D topological insu- 
lators necessarily possess massless Dirac fermion sur- 
face states^, and the fact that they are observed in real 
materials^ have stimulated lots of interest recently. In 
considering the surface electronic structure of such ma- 
terials it is important to remember that the surface of a 
bulk sample is a closed two dimensional manifold. For ex- 
ample, the surface of a slab has the topology of a sphere 
(Fig. (l)a). If we apply an upward magnetic field the net 
flux passing through the surface is zero. The situation 
is shown in Fig. (l)(a) where the magnetic field at the 
top surface points along the surface normal (red), while 
that at the bottom points opposite to the surface normal 
(green). Since there is no magnetic flux piercing the four 
side faces, they are left white. Fig. (l)(b) is topologically 
equivalent to Fig. (l)(a), and the presence of azimuthal 
symmetry makes it a convenient geometry for theoretical 
studies. Because the net magnetic flux is zero there is no 
topologically protected zero energy states^. This raises 
questions about the integrity of the Landau levels and 
the existence of the quantized Hall effect. 

What is the electronic structure of a topologically in- 
sulator with the shape of Fig. (l)(a)? Since the Landau- 
level states on the top and bottom surfaces are connected 
by delocalized states on the side faces, it is not clear what 
the answer should be. It is well known that in an infinite 
plane each Dirac point contributes (n + l/2)e^//i to the 
Hall conductance. What should the Hall conductance of 
a slab be? In addition, because a closed surface has no 
edges, the validity of the conventional edge state picture 
for integer quantum Hall effect^ is questionable even if 
the quantized Hall effect does exist. Finally if, for 5 = 0, 
there is one Dirac point on the surface of a semi-infinite 
sample, how many Dirac points should the surface of a 
slab have? These questions motivate us to look at the 
surface electronic structure of a finite topological insula- 
tor more closely. 

In the following we shall focus on the simplest case 
where the surface Dirac point is derived from of the bulk 
electronic structure near k = (0, 0, 0) (interestingly this is 
the situation of Bi2Se3^). We begin by asking the follow- 
ing question. "If H = Xl/c=i 2^k{—^dk — Ak) describes 
the surface electronic structure of an semi-infinite topo- 
logical insulator what is the surface Hamiltonian for a 
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FIG. 1: (color on-line) The surface of a slab (a) and its topo- 
logical equivalent - the surface of a sphere (b). The red and 
green colors indicate the signs of the magnetic flux being 
opposite. The white region has no magnetic flux piercing. 
(c),(d) the current induced by flux ramping. In part (d) the 
dashed circles are drawn for keeping track of the charge that 
flows pass it. 



finite slab-like sample?" Here k = 1,2 label the two 
orthogonal coordinate frames in the surface; is the 
component of the electromagnetic vector potential; <j are 
the Pauli matrices, and {h/2)n x a, where n is the sur- 
face normal, is the in-plane physical spin operator. For 
simplicity, we have set the Fermi velocity to unity. A 
clarification is in order. For a real material, crystalline 
anisotropy usually results in anisotropy in the low energy 
effective theory. For example the bulk effective Hamilto- 
nian given in Ref.^ clearly exhibits a different Fermi ve- 
locity in the z direction from that in the x, y directions. 
However, this difference can be removed by a rescaling of 
the z coordinate. 

We note that any surface with the topology of a sphere 
has total Gauss curvature 47r hence is necessarily curved. 
It can be shown, by, e.g., projecting the bulk effective 
Hamiltonian in Ref.^ onto surfaces of different orienta- 
tions, that the "spin" a of the surface states is always 
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locked within the surface plane. Due to this spin-space 
locking, the appropriate Hamiltonian is the Dirac Hamil- 
tonian in curved spaces. 

The curvature of the space enters through a non- 
abelian spinor connection^. (Note in our case only the 
space is curved not the time.) In the following we shall 
choose the most convenient coordinate system to write 
down such a Hamiltonain. It has been proven^ that for 
any Riemann surface it is possible to choose a coordinate 
system such that the metric is conformally related to the 
Euclidean metric everywhere, i.e.. 
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Such a coordinate system is called "isothermal" . In terms 
of isothermal coordinates the Dirac Hamiltonian read 



H 



./3(x) 



afe(x) {-idk + irfc(x) - zAfc(x)) . (2) 



fe=l,2 



Here cr/e(x) = e/e(x) • a with ei,2(x) being the local or- 
thonormal coordinate frame so that 63 (x) = ei(x) x 
62 (x) = n is the unit surface normal. The spinor con- 
nection rfc(x) is given by 



and a{u) = 2 tan~''^(e^). As ^ ±00 the point r{u,(j)) 
approaches the south and north pole, respectively. Under 
this coordinate system 



-2f3{z,z) 



(4) 



where z = u -\- and the Gauss curvature is given by 
G{z^z) = Ae^^dzdzP. In the same complex notation the 
magnetic field is given hy B = —2ie^^{dzAz — dzAz). 
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FIG. 2: (color on-line) The eigen spectra of Eq. (7) under 
constant magnetic field - (a), and zero magnetic field - (b). 



rfc(x) = - 



^e/ej5j^(x)e3(x) + ^e3(x) x dkh{^) 
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(Note that due to the choice of isothermal coordinates 
we have used the usual Euclidean space notion which 
does not distinguish the covariant and contravariant in- 
dices.) If we define |z±(x)) as the eigenvectors of 5*3 (x), 
i.e., 5'3(x)|2:±(x)) = ±|^±(x)), and sandwich the Hamil- 
tonian in Eq. (2) between them, the following 2x2 matrix 
operator is obtained^ 



H 



{d-z 





iAz) + ^dzf3 



{dz - iAz) + ldzf3 







(3) 



In the above equation z = x'^ -\- ix^, z = x'^ — ix^, dz = 
l/2{di - 182), d-z = l/2{di + id2), and A, = l/2(Ai - 
iA2)^Az = 1/2(741 +^^2), respectively. In the following 
we shall solve the eigenvalue problem with the magnetic 
field and the configuration space shown in Fig. (l)(b) and 
the Hamiltonian given by Eq. (3). 

First, we choose the isothermal coordinate system 
(ii, (j)) where each point r G 6*^ is parameterized as 
r{u^(j)) = R{sma{u) cos(j)^sma{u)sm(j)^cosa{u)) . Here 
R is the radius of the sphere, — (X)<i^<oo,0<(/)<27r 



If the magnetic field is uniform (hence requires a mag- 
netic monopole to generate it) the field strength is re- 
lated to G hy B = n^G/2 where is the number of 
flux quanta piercing through the sphere. If we pick the 
gauge diAi + 82 A2 = 0, a A{z, z) can be found so that 
Az = -idzK A-z = idzA and B = Ae'^^dzdzA. Gom- 
pare this equation with G = Ae^^dzdzP and B = n0G/2 
suggests that we can set 



(5) 



By substituting Az = —idzA^ Az = idzA and Eq. (5) into 
Eq. (3) we obtain. 



^2 ^ ' 



dz^{l^n^)dzP 







Because of Eq. (4) the Hamiltonian in Eq. (6) is (j)- 
translation invariant, the eigenfunctions are of the form 



r 



i), where m, the azimuthal quantum 



number, is an integer. The equation satisfied by Um and 



Vm IS 



J 



p/3(n) 



E 



(7) 



In the absence of magnetic field (n^ = 0) the eigenspec- trum is shown as Fig. (2)(b). At m = there is one 
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FIG. 3: (color on-line) (a) The eigen spectrum of Eq. (7) 
with magnetic field shown in Fig. (l)(b). The red dashed line 
denotes the Fermi level, (b) \um{u)f -\-\vm{u)f as a function 
of u and m for the upper n = band (the band intersected 
by the Fermi level). Insert of part (b) the The current (J^) 
expectation value associated with the (split) n = Landau 
band. 



pair of eigenvalues with E = 0^ thus there is only one 
Dirac point! For >> 1 the eigen spectrum of Eq. (7) 
consists of energy levels which are independent of m: 
En = ±Y^n(n0 + n)/W' where n > as shown in the 
inset of Fig. (2) (a). This is the relativistic Landau level 
on a sphere^. 

Next, we consider the magnetic field illustrated in 
Fig. (l)(b). Here we choose 

A(^i) = ^(tanh-)/3(^.), (8) 

where w is proportional to the width of the zero-flux band 
in Fig. (l)(a). For the results shown in Fig. (3) w is set 
equal to the magnetic length. Fig. (3) (a) is the eigen 
spectrum. For large m the Landau levels in Fig. (2) (a) 
is recovered. However unlike Fig. (2) (a) where the en- 
ergy levels are non-degenerate at each m, in Fig. (3) (a), 
as m oo, each energy level is two fold-degenerate. 
This degeneracy arises from the Landau states in the 
two hemispheres. In the presence of the Zeeman term 
the Landau level will be split even as m ^ oo. Usually 
that splitting is much smaller than the Landau gap hence 
will not affect any of the discussions below. 

In the following let us focus on the n = Landau level. 
As m decreases the degeneracy are split due to the mix- 
ing with the delocalized states around the equator. Since 
the overlap between the Landau states and the delocal- 
ized states vanishes exponentially as m ^ oo so does the 
splitting. However, strictly speaking so long as m < oo 
the splitting is always non-zero. This is consistent with 
the fact that there is no topologically protected zero 
modes. Spatially the eigenfunctions are localized in the 
zero-flux band around the equator when m is sufficiently 
small (or negative). At large m the eigen functions are 
equal amplitude superpositions of two components local- 
ized in different hemispheres (Fig. (3)(b)). This is due to 



the mixing discussed above. In the inset of Fig. (3)(b) 
we show the expectation value of the current operator, 
^ for large m and develop positive/negative values 
for the upper/lower band as m decreases. If the Fermi 
level is marked by the dashed red line in Fig. (3) (a), it 
will intersect the chiral, current-carrying states. These 
are the "edge" states of our problem. Their dispersion 
are remarkably similar to the valley-split edge states in 
graphene^^. The current in the chiral state is distributed 
across the width of the zero-fiux band. Thus we have 
chiral sheets. 

Now we switch to the existence of quantum Hall ef- 
fect. Let us imagine inserting a thin solenoid as shown in 
Fig. (l)(d) (or (c)) and slowly ramping its flux from to 
00- An adiabatic evolution, m ^ m — 1, will be induced. 
At the end, an extra electron will appear above the Fermi 
level. At the mean time the charges reside in each large m 
states march toward the zero-flux band around the equa- 
tor. Since these states have 50% weight in each hemi- 
sphere (Fig. (l)(a)), each dashed circle in Fig. (l)(d) will 
register an average 1/2 of an electron crossing it. This 
suggests that each hemisphere transfers 1/2 of an electron 
to the equator. The fact that an azimuthal electric field 
induces a quantized charge accumulation in the perpen- 
dicular direction suggests the existence of quantized Hall 
effect. The fact 1/2 + 1/2 = 1 electron is transferred im- 
plies cFxy = /h. If the Fermi level intersects the lower 
n = band, the direction of current will be reversed, 
hence (Jxy = —e^/h. Similar argument carried out for 
higher Landau levels gives cjxy = (2n + l)e^ /h where n 
are integers. These are exactly half of the quantized Hall 
conductances observed in graphene^^. The above result is 
consistent with the statement that dxy = ±(n + l/2)e^//i 
for the two hemispheres. (Here the sign of the Hall con- 
ductance is defined in reference to the surface normal so 
that j = axyfi X E). The fact that the induced EMF 
circulates the surface normal in opposite senses in the 
two hemispheres, and the fact that the Hall conductances 
have opposite sign, imply the induced Hall current will 
flow toward (or away from) the equator in both hemi- 
spheres. Because \(Txy\ = {n-\- l/2)e^//i, the total charge 
that flows to/from the equator is 2 x (n + 1/2) = 2n + 1 
electrons. 

The chiral sheet current discussed above also car- 
ries spin polarization. The effective field theory 
describing the chiral state is given by Hchirai = 
—i J d(j)x^ {^)d(f)X{^)- III the slab geometry (j) is the co- 
ordinate around the four side faces. The spinless fermion 
operator x(0) is related to the original spinful fermion 
operator ^l)^{u, (j)) via x(0) = ! duY,^ ^o,a{u, (l))i^a{u, (j)). 
In the above equation ^o,cr('^^, 0) = Umo{u){(j\zj^{u^(j))) + 
Vmo{u){(7\z-{u,(t))), where {umo{u),Vmo{u)) is the so- 
lution of Eq. (7) for the upper upper n = band 
at m = mo, and \z±{u^(j))) are the eigen spinors of 
(73(1^,0) = 63(1^,0) • a. In Fig. (4)(b) we plot the expec- 
tation value (cr^) = J du 62(1^, (/>) • {^o{u^ (j))\a\£^o{u^ (j))) = 
i J duVmo{u)umo{u) + c.c as a function of mo- Clearly, 
these chiral states also carries a a polarization parallel to 
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FIG. 4: (color on-line) (a) The spectrum of the configuration 
Fig. (l)(b) for w (the width of the zero flux band) equal to the 
magnetic length, (b) The a polarization of the chiral sheet 
states as a function of mo. (c) A two terminal device. Here 
the arrows indicate the direction of charge as well as spin 
currents. (d) The spectrum of the configuration Fig. (l)(b) for 
w equal to ten magnetic length, (e) The set up for measuring 
the quantized Hall effect in the limit thick slab. 



its current direction. However, because the physical spin 
is (h/2)h X (J, the actual spin polarization will be per- 
pendicular to the current direction. This is illustrated 
in Fig. (4)(c), where the direction of the average spin 
polarization and the charge current on the side surfaces 
are indicted by the red and black arrows, respectively. 
Thus in the case of thin slab (i.e. thickness ~ magnetic 
length) a two-terminal measurement (Fig. (4)(c)) will ex- 



hibit both quantized charge conductance and spin Hall 
effect. 

So far we have set the width of the zero-flux band to the 
magnetic length. As the width increases, the edge states 
quickly merge into a continuum, as shown in Fig. (4)(d) 
when w = 10/^- However so long as the Fermi energy 
lies between the n and n + 1 Landau level away from 
the zero-flux band (see the red line in Fig. (4)(d)) the 
previous flux ramping argument still implies the flow of 
2n + 1 electrons into/from the zero-flux band upon the 
adiabatic evolution m m + 1. This suggests that in 
midst of the continuum states, there are a net 2n + 1 
chiral current-carrying edge channels. Due to the anti- 
localization effect^^, both the chiral and the continuum 
states will be delocalized. In this case the two-terminal 
conductance will no-longer be quantized. In order to 
observe the quantized Hall effect one has to set up a 
"Corbino" geometry as shown in Fig. (4)(e). In this case 
the circulating current on the side surface (indicated by 
the red arrow) will be (2n + l)e^ /h times the voltage be- 
tween the electrode and the outer surface of the sample. 

In summary, we solve the curved space Dirac equation 
to study the surface electronic structure of topological 
insulators. For a slab-like sample with magnetic field 
perpendicular to its top and bottom surfaces, there are 
chiral states delocalized on the four side faces. These 
"chiral sheets" carry charge and spin currents. We dis- 
cuss the experimental to measure the quantized charge 
Hall effect {cFxy = (2n + l)e^//i) in the limit of thin and 
thick slabs. Due to the spin-orbit interaction the charge 
current carries spin polarization. In thick slabs, due to 
the anti-localization effect both the chiral and non-chiral 
electronic states will be extended on the side surfaces. 
These states will exhibit unconventional Aharonov-Bohm 
effect in weaker field, when phase coherence is achieved. 
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